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Motivated by the discovery of superconductivity in boron-doped (B-doped) diamond, we 
investigate the localization and superconductivity in heavily doped semiconductors. The com- 
petition between Anderson localization and s-wave superconductivity is investigated from the 
microscopic point of view. The effect of microscopic inhomogeneity and the thermal fluctuation 
in superconductivity are taken into account using the self-consistent 1-loop-order theory with 
respect to superconducting fluctuation. The crossover from superconductivity in the host band 
to that in the impurity band is described on the basis of the disordered three-dimensional attrac- 
tive Hubbard model for binary alloys. We show that superconductor-insulator transition (SIT) 
accompanies the crossover. We point out an enhancement of Cooper pairing in the crossover 
regime. Further localization of the electron wave function gives rise to incoherent Cooper pairs 
and the pseudogap above T c . A global phase diagram is drawn for host band superconductiv- 
ity, impurity band superconductivity, Anderson localization, Fermi liquid state, and pseudogap 
state. A theoretical interpretation is proposed for superconductivity in the doped diamond, 
SiC, and Si. 

KEYWORDS: Anderson localization; superconductor-insulator-transition; localized Cooper pairs; heavily 
doped semiconductors 



1. Introduction 

Superconductivity in the vicinity of quantum phase 
transitions has attracted much interest. Superconductiv- 
ity in high-T c cupratcs, 1 organic superconductors, 2 and 
heavy fermions 3 occurs in proximity to various quantum 
critical points, which has been a central issue in the con- 
densed matter physics. Another classical issue is the SIT, 
which arises from the competition between Anderson lo- 
calization 4 and s-wave superconductivity. Vast theoreti- 
cal studies have been devoted to the SIT triggered by dis- 
orders. 5-16 In this paper, we point out that the recently 
discovered superconductivity in B-doped diamond 17 and 
related materials 18-20 provides a new field of supercon- 
ductivity in the vicinity of Anderson localization. 

Diamond is recognized as a precious gem, but becomes 
a wide-gap semiconductor with light doping with sub- 
stitutional boron acceptors. Superconductivity has been 
discovered by Ekimov et al. by increasing the concen- 
tration of boron acceptors. 17 A transition temperature 
higher than 10 K has been realized in heavily B-doped 
diamond. 21,22 It seems to be surprising that such a high 
transition temperature is induced by low carriers in a 
three-dimensional system since a small density of states 
(DOS) near the band edge is harmful for superconduc- 
tivity. It is highly desired to identify a novel mechanism 
of the enhancement of superconductivity in B-doped di- 
amond. 

Semiconductors have been a central field in exper- 
imental studies of Anderson localization. 23 Low carri- 
ers in semiconductors are significantly affected by ran- 
domness since the electron wave function near the band 
edge is localized by a weak disorder. 24 The effects of 
randomness have been actually revealed by experiments 
on B-doped diamond, which have shown the localiza- 
tion effect, 25 high residual resistivity, 25,26 and absence 



of a well-defined Fermi surface. 27 Superconductivity in 
B-doped diamond is referred as "superconductivity with 
no Fermi surface". 28 However, the localization effect has 
never been investigated in theories of the B-doped dia- 
mond. 29-35 The purpose of this study is to investigate the 
relationship between localization and superconductivity 
in doped semiconductors. 

Two theoretical scenarios have been proposed for 
superconductivity in B-doped diamond. The first one 
is based on the impurity band formed by impurity 
states whose wave functions are localized around boron 
atoms. 29-31 ' 36 The other is based on the host band of 
diamond in which the carriers are provided by substitu- 
tional boron acceptors. 32-35 The former was proposed by 
Baskaran 29 immediately after the discovery of supercon- 
ductivity. The latter has been proposed on the basis of 
first-principles band calculations. 32-35 

Some experimental studies have been carried out to 
determine which scenario well describes the electronic 
structure of B-doped diamond. 25-27,37-39 Although pho- 
toemission spectroscopy 25,27 has never observed any in- 
dication of the impurity band, the presence of the im- 
purity band has been claimed by optical measurement. 38 
We point out here that the two scenarios, namely, super- 
conductivity in the impurity band and that in the host 
band, are not contradicting concepts, and can be con- 
tinuously described using a microscopic model. In this 
study, we examine the crossover from the host band to 
the impurity band on the basis of the disordered attrac- 
tive Hubbard model of binary alloys. We show that the 
crossover is accompanied by the SIT, increases in T c , and 
pseudogap. The crossover in the electric DOS can be 
described by coherent potential approximation (CPA), 
which has been applied to the analysis of B-doped dia- 
mond. 31 However, the spatial inhomogeneity and fluctu- 
ations of superconductivity, which are neglected in the 
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CPA, play an essential role in superconductivity in the 
crossover and impurity band regimes. 
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Table I. Superconducting properties of the B-doped diamond, 
SiC, and Si as well as the Al-doped SiC. The maximum val- 
ues of T c , and H c 2 or H c are shown. More detailed comparisons 
have been given in ref. 40. 

Boron acceptors also induce superconductivity in Si 18 
and SiC. 19 Recently, superconductivity has also been re- 
alized in Al-doped SiC. 20 The superconducting prop- 
erties of these materials are summarized in Table. I. 
More detailed comparisons between B-doped diamond, 
SiC and Si have been summarized by Kriener et al. i0 
The maximum T c values are 0.35 K for B-doped Si and 
1.5 K for B-doped SiC. A more pronounced difference 
from B-doped diamond is the magnetic response in the 
superconducting state. B-doped diamond is a type II su- 
perconductor, while the others are type I superconduc- 
tors. The large Ginzburg-Landau (GL) parameter k ~ 18 
of B-doped diamond 26 should be contrasted to the small 
GL parameter k ~ 0.35 of B-doped SiC. 40 The upper 
critical field of such diamond is 10 3 times higher than 
the critical field of such SiC. These differences should be 
related to electronic properties. B-doped diamond is a 
"bad metal" without any well-defined Fermi surface, 27 
while B-doped SiC is a "good metal" with a small resid- 
ual resistivity 19,40 We will propose here a theoretical 
interpretation for the differences between these super- 
conductors. 

The paper is organized as follows. The disordered 
attractive Hubbard model is introduced in §2. The 
crossover from the impurity band to the host band is 
demonstrated in §3. We examine the Anderson transition 
in §4, and the superconductivity is investigated in §5. We 
discuss B-doped diamond, SiC, and Si in §6. The similar- 
ity between the doped semiconductors and the high-T c 
cuprates is pointed out in §7. These results are summa- 
rized in §8. 

2. Disordered Attractive Hubbard Model 

For a discussion of superconductivity in heavily doped 
semiconductors, we investigate the disordered attractive 
Hubbard model in three dimension: 

H = -t c L c j° + ^2( w i - mK* + U^2rii,im,i, 

<i,j>,<T i i 

(1) 

where rii ia is the number of electrons at the site i with 
the spin a and n, = £ CT nj l<T . The simple cubic lattice is 
assumed and the symbol < i,j > denotes the summation 
over the nearest neighbor sites. We choose the unit of 
energy as t = 1. 

The disorder is described by the random potential Wj. 
We assume a binary alloy in which Wi = at host sites 



(carbon sites in B-doped diamond) and Wi = CA mp at im- 
purity sites (boron sites in B-doped diamond) . We denote 
impurity concentration as nj mp . The chemical potential 
(i is chosen so that the total electron concentration is 
n = 2 — nj m p in accordance with those in B-doped dia- 
mond, SiC, and Si. Note that the presence of a B-H com- 
plex 41 hardly affects the following results, and therefore 
n- lmp should be regarded as the concentration of isolated 
boron atoms. We assume rn mp — 0.02 in §3, §4, and §5.2- 
5.5, and discuss the doping dependence of superconduc- 
tivity in §5.6. 

It is believed that the electron phonon interaction gives 
rise to superconductivity in B-doped diamond. 32-35 We 
take into account the static attractive interaction U < 
to describe the s-wave pairing interaction for simplicity. 

3. Crossover from Host Band to Impurity Band 

The crossover from the host band to the impurity band 
is described by varying the impurity potential [/ imp in the 
model eq. (1). Figure 1 shows the DOS's for various im- 
purity potentials. Here, we assume U = for simplicity. 



(a) U imp =3 (b) U imp =4 




Fig. 1. DOS in the normal state (solid line). We assume U lmp = 
3, 4, 4.6, and 6.1 in (a), (b), (c), and (d), respectively. We assume 
"imp = 0.02 and U = 0. The dashed and dash-dotted lines show 
the partial DOS's in the boron and carbon sites, respectively. 
The numerical calculation is carried out in 21 3 sites. 50 samples 
are taken for determining the random average. 



For a large impurity potential U nup = 6.1, the impurity 
band is clearly separated from the host band, as shown 
in Fig. 1(d). The holes are nearly half-filled in the im- 
purity band. 31 On the other hand, the impurity band is 
implanted into the host band for a small impurity poten- 
tial L7 imp = 3, as shown in Fig. 1(a). Single-particle states 
around the Fermi level are formed by impurity states in 
the former, while those are formed by the host band of di- 
amond in the latter. A crossover between the two regimes 
occurs, and the impurity band merges with the host band 
at around J7 imp = 4 ~ 4.6, as shown in Figs. 1(b) and 
1(c). The acceptor level of the lightly doped boron atoms 
in diamond has been estimated to be 0.37 eV, 42 which 
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corresponds to Ui mp = 4. 6. 31 Thus, B-doped diamond 
seems to be near the crossover regime. 

4. Anderson Localization 

In this section, we investigate Anderson localization in 
the model eq. (1) to clarify the electronic structure in the 
normal state, which is a background of superconductiv- 
ity. We analyze here the inverse participation ratio (IPR) 
and level statistics for this purpose. 

For U — 0, the Hamiltonian eq. (1) is diagonalized on 
the basis of the single-particle states. The IPR is defined 
as 



IPR=<J2\M?1\ 4 >e f , 



(2) 



where ipai'r) is a single-particle wave function with an 
energy e a . The bracket <>e f indicates the average for 
the states around the Fermi energy, \e a — fx\ < e c . We 
choose e c = 0.2 so that the cutoff e c is less than the width 
of the impurity band. The IPR describes the itinerant or 
localized character of the single-particle state. The size 
scaling of the IPR shows IPR oc A^ 1 with N being the 
number of sites for the completely extended state, while 
the IPR is constant for the localized state. We denote 
here N = L 3, with L being the linear size of the three- 
dimensional lattice. The single particle wave function has 
a fractal dimension d* in disordered systems, where IPR 
oc 

N -d'/d_ Note that d 

= 3 in our case. The universal 
distribution of the fractal dimension d* at the mobility 
edge has been shown. 43 The average is estimated to be 
d* = 1.3 for the three-dimensional Anderson transition 
in the orthogonal universality class. 43 
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Fig. 2. Size scaling of the IPR for Ui n 
the bottom to the top. We assume that n- u 



4,4.6,5, and 6.1 from 
: 0.02 and (7 = 0. 
The power against the size N is shown as a visual guide. We 
take into account 220, 150, 100, 70, 50, 40, 30, and 20 samples 
for 13 3 , 15 s , 17 3 , 19 3 , 21 3 , 23 3 , 25 3 , and 27 3 sites, respectively. 



Figure 2 shows the scaling plot of the IPR x TV for a 



fixed impurity concentration m, 



0.02. We observe a 



significant decrease in the fractal dimension d* with in- 
creasing impurity potential U- lmp through the crossover 
from the host band to the impurity band. Since the aver- 
age fractal dimension d* = 1.3 at the mobility edge leads 
to the scaling IPR xN cx TV - 57 , we find Anderson tran- 
sition at around t/i mp = 4.6. Thus, the crossover from 



the host band to the impurity band is accompanied by 
Anderson localization. It is reasonable that the single- 
particle wave function becomes localized with approach 
to the impurity band regime. 

The IPR is also an important index for the supercon- 
ductivity, because the attractive interaction between the 
time-reversed states is proportional to the IPR. Since 
Cooper pairing is determined by this interaction in the 
localized limit, 11 the data in Fig. 2 implies the enhance- 
ment of Cooper pairing with increasing U mip . This is 
qualitatively true, although the parameters assumed here 
are far from the localized limit. The DOS shown in 
Fig. 1 gives another interpretation of the enhancement 
of Cooper pairing across the crossover. The effect of the 
attractive interaction is indicated by the ratio |?7|/W 
for a small C/ imp , while that should be represented by 
l^l/W^imp for a large C/ imp , where W and Wi mp are the 
width of the host band and that of the impurity band, 
respectively. Since W 3> W^imp for a small impurity con- 
centration n- lmp , the interaction is effectively enhanced 
in the impurity band region. The enhanced pairing in- 
teraction leads to the enhancement of superconductivity 
in the crossover regime, as shown in §5. 
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Fig. 3. Results of the level statistics for rii mp = 0.02. We show the 
nearest-level-spacing distribution function P(s) for energy levels 
within \e a — fi\ < e c = 0.2. The crossover from the Wigner- 
Dyson statistics to the Poisson statistics occurs as the impurity 
potential C/i mp increases. The numerical calculation is carried 
out for 21 3 sites and 250 samples. 



Anderson localization in the crossover regime is also 
observed by the analysis of level statistics. 44 Figure 3 
shows the nearest-level-spacing distribution function 
P(s) defined as 



P(s)=<8(\e a 



'-a+l\ 



S) >E F 



(3) 



where e s is the averaged nearest level spacing. It has been 
shown that the distribution obeys the Poisson statistics 
P(s) — exp(— s) in the insulating state, while it obeys the 
Wigner-Dyson statistics P(s) = \ s exp(— 7rs 2 /4) in the 
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metallic state. 44 ' 45 As shown in Fig. 3, the nearest-level- 
spacing distribution function P(s) shows a crossover 
from the Wigner-Dyson statistics to the Poisson statis- 
tics with increasing U lmp across the crossover from the 
host band to the impurity band. Because of the small 
DOS near the band edge, it is difficult to carry out 
finite-size scaling. 44 However, it is clear that the Ander- 
son transition occurs in the crossover regime at around 
Ui mp — 4.6, consistent with the analysis of IPR. 

5. Superconductivity 

5. 1 Formulation 

Next, we discuss the superconductivity. We adopt real- 
space self-consistent T-matrix approximation (RSTA) to 
analyze the model eq. (1). The RSTA has been formu- 
lated for disordered <i-wave superconductivity in high-T c 
cuprates. 46 The thermal fluctuation of superconductivity 
is taken into account in the self-consistent 1-loop order. 
Because the calculation is carried out in real space, ran- 
domness is accurately taken into account. 

Most microscopic studies of superconductivity in ran- 
dom systems have been based on mean field approxi- 
mation. 10, 47-49 However, the fluctuation plays an impor- 
tant role in strongly disordered systems because the fluc- 
tuation is enhanced by microscopic inhomogeneity. 10,46 
RSTA can be carried out in more than 10 3 sites, which 
is needed for the analysis of three-dimensional models. 
We carry out a numerical calculation in a 11 x 11 x 11 
lattice, unless otherwise mentioned. This lattice size is 
much larger than the limit of quantum 12 and classical 13 
Monte Carlo simulations. 

The formulation of RSTA is described as follows. We 
diagonalize the non-interacting Hamiltonian 

H = -t J2 clc^ + Y^iWi-^rii, (4) 

<2, i 7>,(T i 

and obtain the undressed Green function as 

G (r,P, w „) = V M^— Kir'), (5) 

a 

where uj n = (2n + 1)ttT is the Matsubara frequency. 
The Dyson equation is expressed as 

G(r, r', u n ) = G (r, r', u n ) + 

2J G (f,fi,uj n )'E(fi,f2,u> n )G(r2,r',ui n ), (6) 

where £(ri,r2,w n ) is the self-energy represented in real 
space. We estimate self-energy using the self-consistent 
1-loop-order approximation, whose results in the clean 
limit have been summarized in ref. 50. We adopt quasi- 
static approximation 51-55 as in ref. 46, where quantum 
fluctuation is ignored, but thermal fluctuation is ap- 
propriately taken into account. We have quantitatively 
shown the validity of quasi-static approximation as well 
as self-consistent 1-loop-order approximation for super- 
conducting fluctuation. 55 Note that the characteristic dy- 
namical property of superconducting fluctuation leads to 
the validity of quasi-static approximation. 55 
Self-energy is given as 

S(r, ?, u n ) = -TU 2 T(f, ?)G{?, r, (7) 



The propagator of superconducting fluctuation, namely, 
the T-matrix, is obtained using 

T(r, ?) = T (f, ?) - UT (f, r x )T{r x ,?), (8) 
fi 

Toft?) = T^G(f,r>„)G(r>Vwn)^W 2 , (9) 

n 

where we take into account the cutoff frequency as 
(j>(uj) = [1 + exp{a(\u\/uj c - l)}] -1 with a = 20. We 
assume the cutoff frequency lu c = 0.2 in accordance with 
the phonon frequency w p h ~ 150 meV in B-doped dia- 
mond. 39 The T-matrix diverges at the superconducting 
transition temperature. This criterion gives rise to the 
transition temperature of the BCS theory, namely, the 
mean field theory with respect to U, when we neglect 
self-energy. 

Since the random potential is included in the unper- 
turbed Hamiltonian Hq, randomness is accurately taken 
into account in RSTA. The spatial dependences of super- 
conducting fluctuation and single-particle states are self- 
consistently determined. We find that perturbative ap- 
proximations for randomness, such as the Born approx- 
imation, 56 CPA, 31 and T-matrix approximation, 57-59 
break down for a large impurity potential L7 imp > 4. 
This is because the quasiparticle states around the Fermi 
level mainly consist of impurity states whose wave func- 
tion is localized around impurity sites. Impurity potential 
should be taken into account in zeroth-order approxima- 
tion when the impurity band plays an essential role in 
superconductivity. 

5.2 Microscopic inhomogeneity 

First, we investigate the spatial inhomogeneity of su- 
perconducting susceptibility defined as 

Xsc(r1 = ^T(r,r0, (10) 

r' 

which is spatially uniform in the clean limit. A typical 
spatial dependence of Xacif) is shown in Fig. 4 for var- 
ious L7i mp values. Since it is difficult to show the spa- 
tial dependence in the entire three-dimensional lattice, 
r = (x,y,z), we show the spatial dependence on the 
plane z = 1. 

(a)Uim P =3 (b)Uim P =4.6 (c)Uim P =6.1 




Fig. 4. Typical spatial dependences of the superconducting sus- 
ceptibility Xsc{r) for (a) l/i mp = 3, (b) Ci mp = 4.6, and (c) 
U inl p = 6.1. We assume n imp = 0.02, T = 0.002, and U = -1. 
We show the results on the plane at z = 1 in the three- 
dimensional lattice. 
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In Figs. 4(a)-4(c), the superconducting correlation de- 
velops in the "dirty region" that includes many impuri- 
ties, in contrast to the conventional spatial dependence, 
in which superconductivity is suppressed by the impuri- 
ties. 10,46-49 This is because the wave functions of quasi- 
particle states around the Fermi level have a large weight 
around impurity sites, and therefore, impurity sites play 
more important roles in superconductivity than host sites 
regardless of Ui mp > 0. Thus, a unique spatial structure 
of superconductivity is realized in heavily doped semi- 
conductors. The s-wave symmetry of superconductivity 
is essential for this result, as discussed in §7. 

The spatial inhomogeneity of the superconducting cor- 
relation, namely, the mesoscopic fluctuation 60 increases 
with increasing Ui mp . Superconductivity is spatially ho- 
mogeneous in the host band regime U[ mp < 4, but sig- 
nificantly inhomogencous in the impurity band regime 
t^imp > 5, as shown in Fig. 4. The microscopic inho- 
mogeneity of the superconducting correlation can be re- 
garded as the localization of Cooper pairs. 11 Figure 4 in- 
dicates that the Anderson localization of single-particle 
states is accompanied by the localization of Cooper pairs. 
A similar microscopic inhomogeneity has been discussed 
for high-T c cuprates from both experimental 61 and the- 
oretical points of view. 13,46-49 ' 62 ' 63 

5.3 Superconducting fluctuation and pseudogap 

The spatial inhomogeneity of the order parameter gen- 
erally disturbs the long-range coherence of superconduc- 
tivity. Then, the short-range correlation induces a pseu- 
dogap in the excitation spectrum. In other words, the 
fluctuation is enhanced by the localization of Cooper 
pairs, and therefore the precursor of superconductivity 
appears above T c . 

Figure 5 shows the DOS above T c , which is obtained 
using 

p(w) =< ~ yimG R (r,f,io) > av , (11) 



(a) Uim P =4, T=0.0007 



(b) Uim P =4.6, T=0.002 



ttN ^ 



where <> av denotes the random average. We take the 
random average for 10 samples of the impurity distribu- 
tion. We find that the DOS is self-averaged, and therefore 
the sample dependence is weak for the spatially averaged 
DOS. We see that the effect of superconducting fluctua- 
tion is negligible in the host band regime ((7i mp = 4), 
while the pseudogap is induced by a short-range su- 
perconducting fluctuation in the crossover and impurity 
regimes (t/i mp = 4.6,5, and 6.1). It is shown that the 
pseudogap increases with increasing impurity potential, 
and therefore the DOS at the Fermi energy is signifi- 
cantly suppressed in the impurity band regime above T c . 

We explain here the mechanism of pseudogap in detail. 
(1) The small width of the impurity band leads to a small 
Fermi energy, and therefore a small coherence length. A 
pseudogap generally appears in small-coherence-lcngth 
superconductors, such as high-T c cuprates and organic 
superconductors. 50 (2) The fractal structure of the 
single-particle wave function enhances the pseudogap, as 
discussed by Feigel'man et al. 11 (3) The localization of 
Cooper pairs, namely, the microscopic inhomogeneity of 
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(c) Uim P =5, T=0.002 





Fig. 5. DOS's for various impurity potentials U- lulp (solid lines). 
We choose the temperature close to the critical point of the su- 



perconductivity, T = 0.0007 for U- n 



4, and T = 0.002 for 



Uimp = 4.6,5, and 6.1. We take the random average over 10 
samples of the impurity distribution. We assume n lmp = 0.02 
and U = — 1. The bare DOS at U = is shown for a comparison 
(dashed lines). 



superconductivity enhances the thermal fluctuation and 
also the pseudogap. 46 (4) The effective pairing interac- 
tion for the s-wave superconductivity is enhanced by the 
localization of single particle-states, as discussed in §4. 
The transition temperature in the mean field theory T C MF 
with respect to U, where the fluctuation is neglected, in- 
creases with increasing impurity potential Ui mp . These 
features give rise to the pseudogap in the crossover and 
impurity band regimes. 

The pseudogap disappears in the host band regime 
because of the wide host band of diamond, the ex- 
tended wave function of quasiparticle states, and the spa- 
tially homogeneous superconducting correlation. Thus, 
the crossover from host band superconductivity to im- 
purity band superconductivity is accompanied by the 
crossover from the conventional Fermi liquid state to the 
pseudogap state. 

A large attractive interaction U = —9.2 gives rise to 
the transition temperature T c ~ 0.002 in the clean limit 
(EAmp = 0). We confirmed that the superconducting fluc- 
tuation is negligible for these parameters. The thermal 
fluctuation at T = 0.002 hardly affects the DOS in the 
absence of disorders. This means that randomness plays 
an essential role in the pseudogap in the crossover and 
impurity band regimes. 

5.4 Superconducting susceptibility 

We investigate here the superconducting susceptibility 
averaged over real space as well as the randomness 



Xs 



< -^^2Xsc(r) > av • 



(12) 



This susceptibility diverges at the critical point of super- 
conductivity. Although the superconducting susceptibil- 
ity is always finite in our calculation due to the finite size 
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effect, the magnitude of Xsc indicates the closeness to the 
critical point. 



O 
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2% doping 
4% doping 




6.5 



Fig. 6. Superconducting susceptibility Xsc- We show the t/j mp de- 
pendences for ni mp = 0.02 (circles) and ni mp = 0.04 (squares), 
respectively. We assume T = 0.002 and U = —I. 



Figure 6 shows the C/; mp dependences of Xsc for n imp = 
0.02 and nj mp = 0.04. We see the dome shape of super- 
conducting susceptibility against the impurity potential 
for a fixed impurity concentration. The peak appears in 
the crossover regime. These results indicate that the op- 
timum T c is obtained in the crossover regime. 
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Fig. 7. Superconducting susceptibility Xsc estimated for the su- 
percell (diamonds). We put the impurities at r = (31 + 1, 3m + 
1, 3n + 1) with /, m, and n being integers. The impurity concen- 
tration is nj mp = 0.0203 for the 11 X 11 X 11 lattice. The Xsc in 
the random system with nj mp = 0.02 is shown for comparison 
(circles). We assume T = 0.002 and U = -1. 



The randomness is essential for the dope shape of 
the superconducting susceptibility in Fig. 6. To illumi- 
nate this point, we show the superconducting suscepti- 
bility x sc in the supercell, in which impurities are set at 
r = (3/ + 1, 3m + 1, 3n + 1). Figure 7 shows the compar- 
ison between the supercell and the random system. We 
see that the susceptibility x sc monotonically increases 
in the supercell with increasing £/j mp , in contrast to the 
random system. Thus, the suppression of the supercon- 
ductivity in the impurity band regime rather arises from 
randomness. 

The enhancement of superconductivity in the supercell 
has been pointed out by Shirakawa et al. 31 on the basis 



of the CPA. Our results are consistent with their claim. 
However, superconductivity is enhanced in the supercell 
in our calculation because fluctuation, which is not taken 
into account in the CPA, is suppressed. 

We have confirmed that the transition temperature in 
the mean field theory, T C MF , monotonically increases with 
increasing impurity potential C/i mp owing to the enhance- 
ment of pairing interaction, as discussed in §4. On the 
other hand, superconducting susceptibility is suppressed 
by mesoscopic and thermal fluctuations in the crossover 
and impurity band regimes. Thus, the competition be- 
tween the enhancement of Cooper pairing and that of su- 
perconducting fluctuations gives rise to the dome shape 
of superconducting susceptibility in Fig. 6. 

The long-range coherence of superconductivity is 
hardly realized in the impurity band regime because of 
the localization of Cooper pairs (see §5.2). The suscep- 
tibility defined in eqs. (10) and (12) takes into account 
both the short-range and long-range correlation of su- 
perconductivity. The weight of the short-range and long- 
range correlation can be clarified by analyzing the aver- 
aged correlation function 



-YT(P 

n 4^ 



+ r,r') > a 



(13) 



We have confirmed that the correlation function in- 
creases for a small \r\ in the impurity band regime, but 
then decreases for a large \r\ with increasing Ui mp . This 
means that the short-range correlation develops, but the 
long-range correlation is suppressed in the impurity band 
regime, as we have mentioned above. 

Figure 6 shows that the peak of Xsc shifts to the im- 
purity band regime with increasing impurity concentra- 
tion rii mp . This is because Anderson localization is sup- 
pressed by increasing the concentration of carriers, n imp , 
and then superconducting fluctuation is also suppressed. 
The pseudogap is suppressed by increasing n- lmpj while 
the T c is increased, as will be shown in Figs. 9-11. 

5.5 Phase diagram 

We illustrate here a schematic figure of the phase di- 
agram against the impurity potential U{ mp and the tem- 
perature T (Fig. 8). 

As we have discussed in §4, the pairing interaction is 
effectively enhanced by the localization of single-particle 
wave functions with increasing U lmp . Therefore, the tran- 
sition temperature T C MF in the mean field theory and the 
magnitude of superconducting gap A at T — monoton- 
ically increase with increasing impurity potential U- lmp 
(dashed line in Fig. 8). We have confirmed this U lmp de- 
pendence of the superconducting gap in the mean field 
theory at T = 0. 

The Anderson localization of single-particle wave func- 
tions leads to microscopic inhomogeneity, which en- 
hances the thermal fluctuation in superconductivity. 
Then, the long-range coherence is disturbed, as shown in 
§5.4. The competition between the enhancement of the 
effective pairing interaction and that of superconducting 
fluctuations leads to a dome shape of the superconduct- 
ing state with a peak T c in the crossover regime (solid 
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Fig. 8. Schematic figure of the phase diagram against the impu- 
rity potential f/i mp and the temperature T. The crossover from 
the host band to the impurity band is described. The solid line 
shows the transition temperature T c of superconductivity. The 
dashed line shows the transition temperature T ( ! vIF in the mean 
field theory. The shaded region indicates the pseudogap state in- 
duced by the incoherent Cooper pairs. The upward arrow (C/mi) 
shows the virtual quantum metal-insulator transition in the ab- 
sence of superconductivity. The downward arrows indicate our 
interpretation for B-doped diamond, SiC, and Si (see §6). 



line in Fig. 8). 

The short-range correlation of superconductivity gives 
rise to the resonance scattering between incoherent 
Cooper pairs and quasiparticles, which leads to a pseu- 
dogap in the shaded region in Fig. 8. In other words, the 
superconducting gap obtained in the mean field theory 
is changed to a pseudogap by the fluctuations. 

The pseudogap state would be an insulating state 
in the highly disordered regime because both electrons 
and Cooper pairs are localized. This state is regarded 
as the "hard gap insulator" proposed by Feigel'man et 
al., which seems to be realized in the disordered thin 
films. 64-75 Thus, the insulating state near the SIT is ex- 
pected to be the localized Cooper pairing state, but not 
the Coulomb glass state. 76 

5.6 Doping dependence 

We investigate here the doping dependence. Impurity 
level depends on the impurity concentration n; mp in the 
heavily doped semiconductors, but we show the results 
for a fixed l/; mp for simplicity. 

Figure 9 shows the doping dependence of the Fourier- 
transformed correlation function 

T*(q) = J T av (r)e^dq. (14) 

The spatially averaged superconducting susceptibility 
Xsc is related to this correlation function as Xsc = T£ V (Q). 
We see that the correlation function increases at q = 
with increasing doping concentration; it hardly changes 
at q 7^ 0. This means that the long-range correlation of 
superconductivity develops in the heavily doped regime. 
This is consistent with the experimental observation that 
the T c of B-doped diamond increases with increasing con- 
centration of isolated boron acceptors. 21,22 ' 41 



6r 




Fig. 9. Fourier-transformed correlation function T^ v (q) along q = 
qx for n imp = 0.01, 0.02, 0.04, 0.08, and 0.125 from the bottom 
to the top. We assume U irnp = 4.6, U = -1, and T = 0.002. 



It has been pointed out that the concentration of car- 
riers is different from that of boron atoms owing to the 
presence of the B-H complex and/or the B-B dimer. 41,77 
We have confirmed that the B-H complex as well as the 
B-B dimer neither enhance nor suppress superconductiv- 
ity. Therefore, the doping concentration n- lmp discussed 
in this paper should be regarded as the concentration of 
isolated boron acceptors, consistent with the experimen- 
tal observation. 41 




CO 

Fig. 10. DOS's for n imp = 0.01, 0.02, 0.04, 0.08, and 0.125 from 
the bottom to the top. We assume that C/i mp = 4.6, U = — 1, 
and T = 0.002. The vertical dashed line is drawn to emphasize 
the doping dependence of the gap in DOS. 



Figure 10 shows the doping dependence of the DOS 
in the crossover regime (Ui mp — 4.6). We see that the 
pseudogap is suppressed by increasing the doping con- 
centration nimp, although the superconducting correla- 
tion is enhanced, as shown in Fig. 9. The single-particle 
states around the Fermi level become itinerant with in- 
creasing n lmp owing to the percolation of impurity states. 
Then, Anderson localization is suppressed, and therefore 
superconducting fluctuation is suppressed. 

The schematic figure of the doping dependence in the 
crossover regime is shown in Fig. 11(a). We would like to 
stress again that the transition temperature T c of su- 
perconductivity increases, but the thermal fluctuation 
is suppressed by acceptor doping. The pseudogap and 
superconducting gap decrease in contrast to T c with 
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(a) Crossover regime (b) Host band regime 




imp 1 'imp 



Fig. 11. Schematic figure of the doping dependence in (a) the 
crossover regime and (b) the host band regime, respectively. We 
show the magnitude of the superconducting gap (or pscudogap) 
divided by 7 = 1.73 for a comparison with T c . 



increasing doping concentration. These features in the 
crossover regime should be contrasted to those in the host 
band regime where the thermal fluctuation is negligible 
except at a very low concentration of carriers. A conven- 
tional doping dependence is expected in the host band 
regime; superconducting gap increases together with T c , 
as shown in Fig. 11(b). 

We have investigated the effect of doping compensa- 
tion due to doped donors and carrier increase due to the 
additional acceptors. We found that superconductivity 
is significantly suppressed. 78 In other words, high T c of 
superconductivity is expected in uncompensated semi- 
conductors with large carrier concentration. 

6. Interpretation for B-doped Diamond, SiC and 
Si 

We discuss here B-doped diamond, SiC, and Si on the 
basis of the results in §5. Since no clear impurity band 
is observed in the angle-resolved photoemission spec- 
troscopy (ARPES), 27 B-doped diamond does not seem 
to be in the impurity band regime. On the other hand, 
the high transition temperature T c > 10 K, 21 ' 22 high up- 
per critical field H c2 > 10 T, 21 ' 22,26 and large Ginzburg- 
Landau (GL) parameter k ~ 18 26 imply that B-doped 
diamond is close to or in the crossover regime. In general, 
the localization of single-particle wave functions leads to 
a small coherence length and a low superfluid density. 
Therefore, a high upper critical field H C 2 and a large GL 
parameter k are expected in the crossover regime. The lo- 
calization effect revealed by electric resistivity measure- 
ments 21 ' 25 also indicates that B-doped diamond is close 
to the crossover regime. 

These observations should be contrasted to those of 
B-doped SiC. B-doped SiC is a type I superconductor 
whose critical field H c is 10 -3 times less than the upper 
critical field of B-doped diamond. 19 Its GL parameter 
K ~ 0.35 19,40 is much less than that of B-doped dia- 
mond. A small residual resistivity implies that the effect 
of Anderson localization is negligible for superconductiv- 
ity in B-doped SiC. These marked differences between di- 
amond and SiC indicate that B-doped SiC is in the host 
band regime far from the crossover regime. This is con- 
sistent with the acceptor level of boron atoms in the SiC 
~ 0.3 eV, 20 ' 79 which is smaller than that in diamond. It 
is expected that the B-doped Si is far from the crossover 
regime because its acceptor level of ~ 0.045 eV 20 ' 80 is 
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the smallest among those of diamond, SiC, and Si. We 
have illustrated these interpretations in Fig. 8. 

Before closing the section, we will comment on the dif- 
ference between B-doped SiC and Al-doped SiC. Both 
compounds have similar T c values of ~ 1.5 K. However, 
Al-doped SiC seems to be a type II superconductor with 
k ~ 1.8 in contrast to B-doped SiC being a type I su- 
perconductor. Since the acceptor level of Al-doped SiC 
of ~ 0.2 eV is smaller than that of B-doped SiC, the 
difference in the GL parameter k seems to be incom- 
patible with the interpretation based on Fig. 8. How- 
ever, the difference between B-doped SiC and Al-doped 
SiC may be understood in the following way. We expect 
that silicon sites play more important roles in the su- 
perconductivity than carbon sites, because the DOS at 
the Fermi level mainly arises from silicon sites. 20 It has 
been pointed out that boron atoms mainly substitute for 
carbon atoms. 20 If aluminium atoms substitute for sil- 
icon atoms rather than for carbon atoms, the random 
potential due to substitutional aluminum acceptors will 
affect superconductivity more significantly than that due 
to boron acceptors. Then, type II superconductivity can 
occur in Al-doped SiC. Further studies are desired to 
examine this possibility. 

7. Similarity to High-T c Cuprates 

We point out the close relationship between doped 
semiconductors and high-T c cuprates. We see the similar- 
ity between Fig. 8 and the typical phase diagram of high- 
Tc cuprates. s-Wave superconductivity and Anderson lo- 
calization in doped semiconductors correspond to d-wave 
superconductivity and the antiferromagnetic Mott insu- 
lator in high-T c cuprates, respectively. The host band 
and impurity band regimes in the former correspond to 
the overdoped and underdoped regimes in the latter, re- 
spectively. 

We have claimed that the pseudogap in the under- 
doped high-T c cuprates is a precursor of superconduc- 
tivity. 50,81 ' 82 The pseudogap due to the same mecha- 
nism has been discussed for the doped semiconductors in 
this paper. The same mechanism of the pseudogap has 
also been proposed to high-density quark matters. 83 Note 
that the pseudogap in these systems does not arise from 
the shift in chemical potential, which describes the BCS- 
BEC crossover in zeroth-order approximation. 84-86 In 
general, the shift in chemical potential is quantitatively 
important in low-density systems, such as trapped cold 
fermion gases, 87 while higher-order corrections, namely, 
resonance scattering between incoherent Cooper pairs 
and quasiparticles, play an essential role in high-density 
systems. 50 Although a doped semiconductor seems to be 
a low-density system, its density is effectively high in 
the impurity band regime because the impurity band is 
nearly half-filled. Actually, the shift in chemical potential 
is negligible in our calculation. The BCS-BEC crossover 
described by the shift in chemical potential was first pro- 
posed for superconductivity in semiconductors. 8 How- 
ever, this is not the case in our calculation. 

Many lines of experimental evidence have been 
obtained for microscopic inhomogeneity in high-T c 
cuprates. 61 According to the results of RSTA, the mi- 
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croscopic inhomogcncity generally appears in small- 
coherence-length non-s-wave superconductors. 46 In con- 
trast to that of the s-wave superconductor discussed in 
this paper, the microscopic inhomogcncity of non-s-wave 
superconductors is not necessarily accompanied by the 
localization of quasiparticles. A weak disorder leads to 
microscopic inhomogeneity in the non-s-wave supercon- 
ductor with a small coherence length. On the other hand, 
the s-wave superconductivity is robust against nonmag- 
netic disorders in accordance with the Anderson theo- 
rem 88 until Anderson localization occurs. 

We note that the microscopic inhomogeneity of high-T c 
cuprates is enhanced by the multi-criticality of antifcr- 
romagnetism and ci- wave superconductivity. 63 An anti- 
ferromagnetic correlation develops in the "dirty region" , 
while a local superconductivity occurs in the "clean re- 
gion". It is expected that the insulating state near the 
SIT in high-T c cuprates is neither the simple Mott in- 
sulator nor the Anderson insulator in analogy with the 
doped semiconductors discussed in §5.5. It is expected 
that the localized Cooper pairing state is realized in a 
variety of systems near the SIT. 

8. Summary and Discussion 

We have investigated localization and superconductiv- 
ity in doped semiconductors on the basis of the disor- 
dered attractive Hubbard model in three dimensions. We 
focused on the crossover from the superconductivity in 
the host band to that in the impurity band. Anderson 
localization is accompanied by the crossover. We found 
that the effective interaction for the paring is enhanced 
in the crossover regime. The superconducting correlation 
is enhanced in the crossover regime, but suppressed by 
the mesoscopic and thermal fluctuations in superconduc- 
tivity with approach to the impurity band regime. The 
dome shape of the superconducting state is expected in 
the phase diagram against the impurity potential U lmp 
and the temperature T. The long-range coherence is de- 
stroyed in the impurity band regime, and then the short- 
range correlation leads to a pseudogap. The insulating 
state due to the localization of Cooper pairs is expected 
to be realized in the highly disordered regime. 

We proposed that the marked differences between B- 
doped diamond, SiC, and Si are understood by taking 
into account the impurity potential of these compounds, 
which are determined by the acceptor level in the lightly 
doped region. B-doped diamond is close to or in the 
crossover regime, while the others are in the host band 
regime far away from the crossover regime. The exper- 
imental results on transition temperature, critical mag- 
netic field, residual resistivity, and the feature of the SIT 
are consistent with our interpretation. 

Two mechanisms have been proposed for the SIT in 
s-wave superconductors. One is the so-called fermionic 
effect, which arises from the long-range Coulomb inter- 
action. 5 ~ 7 ' 9 The other is the bosonic effect, namely, the 
fluctuations in superconductivity. 8, 10, n ' 13-16 Although 
the SIT in some bulk materials can be understood by tak- 
ing the fermionic effect into account, 24 the bosonic effect 
has been indicated by the recent experiments on disor- 
dered thin films, such as Bi, InO^, and TiN films. 16 ' 66 - 75 



We have pointed out that a significant bosonic effect ap- 
pears in a doped wide-gap semiconductor, which may be 
realized in B-doped diamond. 

In summary, we have investigated superconductivity 
and localization with focus on the crossover from the 
host band regime to the impurity band regime in doped 
semiconductors. The superconducting correlation is en- 
hanced in the crossover regime where the electrons have 
a dual nature, namely, the itinerant and localized char- 
acters. Anderson localization in the crossover regime is 
accompanied by the mesoscopic and thermal supercon- 
ducting fluctuations in superconductivity. The similarity 
of doped semiconductors to high-T c cuprates and disor- 
dered thin films has been discussed. 
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